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Abstract 

Consider the p-system describing the subsonic flow of a fluid in a pipe 
with section a — a(x). We prove that the resulting Cauchy problem 
generates a Lipschitz semigroup, provided the total variation of the 
initial datum and the oscillation of a are small. An explicit estimate 
on the bound of the total variation of a is provided, showing that 
at lower fluid speeds, higher total variations of a are acceptable. An 
example shows that the bound on TV (a) is mandatory, for otherwise 
the total variation of the solution may grow arbitrarily. 

2000 Mathematics Subject Classification: 35L65, 76N10. 
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1 Introduction 

Consider a gas pipe with smoothly varying section. In the isentropic or 
isothermal approximation, the dynamics of the fluid in the pipe is described 
by the following system of Euler equations: 



d t (ap) + d x (aq) = 

,,2 

d t (aq) +d x 



a I — +PO) 



P (p) d x a, 



(1.1) 



where, as usual, p is the fluid density, q is the linear momentum density, 
p = p(p) is the pressure and a = a(x) is cross-sectional area of the tube. 
We provide a basic well posedness result for (|l.lj) . under the assumptions 
that the initial data is subsonic, has sufficiently small total variation and 
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the oscillation in the pipe section a = a{x) is also small. We provide an 
explicit bound on the total variation of a. As it is physically reasonable, as 
the fluid speed increases this bound decreases and vanishes at sonic speed, 
see dHU). 

As a tool in the study of (|l.lj> we use the system recently proposed for 
the case of a sharp discontinuous change in the pipe's section between the 
values a~ and a + , see [31 El E] ■ This description is based on the p-system 



d t p + d x q = 



(1.2) 



equipped with a coupling condition at the junction of the form 

*(o-,(p, g )(t J 0-);o + ,(p,g)(t,0+)) =0 (1.3) 

whose role is essentially that of selecting stationary solutions. 

Remark that the introduction of condition (|1.3p is necessary as soon as 
the section of the pipe is not smooth. The literature offers different choices 
for this condition, see [3j[71[8]. The construction below does not require any 
specific choice of ^ in (jl.3p . but applies to all conditions satisfying minimal 
physically reasonable requirements, see (S0)-(X12) 

On the contrary, if a £ W ' the product in the right hand side of the 
second equation in (jl.ip is well defined and system (jl.ip is equivalent to the 
2x2 system of conservation laws 

dtp + d x q = — | d x a 

d tq + d x {^+p(p))=-^d x a. (L4) 

Systems of this type were considered, for instance, in [Til EE US OS 
HSimj. In this case the stationary solutions to (jl.ip are characterized as 
solutions to 

d x (a(x)q) = (g xq = _l dxa 

d x (a(x) (£+p(p))) =p{p)d x a OT \d s (£+p{p))=-£d x a, (L5) 

see Lemma 12.61 for a proof of the equivalence between (jl.4p and (jl.ip . 

Thus, the case of a smooth a induces a unique choice for condition (II. 3D . 
see (j2.3p and (j2.19p . Even with this choice, in the case of the isothermal 
pressure law p{p) = c 2 p, we show below that a shock entering a pipe can 
have its strength arbitrarily magnified, provided the total variation of the 
pipe's section is sufficiently high and the fluid speed is sufficiently near to 
the sound speed, see Section 12.21 Recall, from the physical point of view, 
that the present situation neglects friction, viscosity and the conservation of 
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energy. Moreover, this example shows the necessity of a bound on the total 
variation of the pipe section in any well posedness theorem for (jl.ip . 

The next section is divided into three parts, the former one deals with a 
pipe with a single junction, the second with a pipe with a piecewise constant 
section and the latter with a pipe having a W 1 ' 1 section. All proofs are 
gathered in Section [3l 



2 Notation and Main Results 



± 



Throughout this paper, u denotes the pair (p, q) so that, for instance, u 
(p^^q^ 1 ), u = (p,q), ■ ■ .. Correspondingly, we denote by f(u) = [q,P(p,q)^j 



the flow in (11.20 . Introduce also the notation R + = [0, +oo[, whereas M + = 
]0,+oo[. Besides, we let a(x±) = lim^ x ± a(£). Below, B(u;6) denotes the 
open ball centered in u with radius 5. 

The pressure law p is assumed to satisfy the following requirement: 



(P) peC 5 



") is such that for all p > 0, p'(p) > and p"(p) > 0. 



The classical example is the 7- law, where p(p) = k p 1 , for a suitable 7 > 1. 

Recall the expressions of the eigenvalues Ai 5 2 and eigenvectors r\p °f the 
p-system, with c denoting the sound speed, 



Ai(u) = ;-c( P ) 

n(u) = 



-l 
Ai(«) 



c(p) = Vp'(p) > A 2<» 
, r 2 {u) ■- 



p 
1 

A 2 (n) 



+ c(p) 



The subsonic region is given by 

A = ju G R+ x R: Ai(u) < < A 2 (n)} 
For later use, we recall the quantities 



(2.1) 



(2.2) 



flow of the linear momentum: P{u) 
total energy density: E(u) 
flow of the total energy density: F(u) 



— +P{p) 
P 



■ p(r) 



dr 



l.(E(u)+p(p)) 



where p* > is a suitable fixed constant. As it is well known, see |10|. for- 
mula (3.3.21)], the pair (E,F) plays the role of the (mathematical) entropy 
- entropy flux pair. 
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2.1 A Pipe with a Single Junction 



This paragraph is devoted to fjl.2[> — fjl.3[> . Fix the section a > A, with A > 
and the state u G A$. 

First, introduce a function E = E(a~,a + ;u~) that describes the effects 
of the junction when the section changes from a~ to a + and the state to the 
left of the junction is u~ . We specify the choice of (|1.3p writing 



arq^ — a q 
a+P(u + ) - a~P(u~) 



S(a",a + ;u~). (2.3) 



We pose the following assumptions on E: 
(SO) E G C 1 ([a — A, a + A] x B (u; 5) ; R 2 ). 

(El) E(a, a; u _ ) = for all a G [a — A, a + A] and all u~ G -B(it; 5). 

Condition (SO) is a natural regularity condition. Condition (SI) is aimed 
to comprehend the standard "no junction" situation: if a~ = a + , then the 
junction has no effects and E vanishes. 

Conditions (SO)-(Sl) ensure the existence of stationary solutions to 
problem <[0]> -(j01l. 

Lemma 2.1 Let (^O)-(Til) hold. Then, for any a G M + , u G Aq, there 
exists a positive 5 and a Lipschitz map 



T: ] a - 5, a + 5 [ x ] a - 5, a + 5 [ x B (u; 5) -> A 



(2.4) 



such that 



^f(a~ ,u~; a + , u + ) 
a~ G a — 5, a + 6 
a + G a — 5, a + 6 
u~ ,u + G B Ou; S) 



T(a , a + ;u ) . 



In particular, T(a, a, u) = u. We may now state a final requirement on S: 

(S2) £(a-,a°;-u-) + E (a , a + ; T(a _ , a ; u - )) = S(a-,a + ;u") . 

With T as in Lemma [2. 11 Alternatively, by (|2.3|) . the above condition (S2) 
can be restated as 



fy(a~,u~;a°,u ) 
y(a°,u°;a + ,u + ) 



*(a ,u ;a + ,u + ) = 0. 



Condition (S2) says that if the two Riemann problems with initial states 
(a - , u~), (a , u°) and (a , u ), (a + , u + ) both yield the stationary solution, 
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then also the Riemann problem with initial state (a~,u~) and (a + ,u + ) is 
solved by the stationary solution. 

Remark that the "natural" choice (|2.19p implied by a smooth section 
satisfies (SO), (SI) and (S2). 

Denote now by u a map satisfying 

v . JiTifx<0 * (a-,u~;a + ,u + ) = 0, . . 

u(x) = _ with ^ i . ; (2.5) 

v ' [n + ifx>0 u~,u + £A . v ' 

The existence of such a map follows from Lemma l2.ll Recall first the defi- 
nition of weak ^-solution, see Definition 2.1] and [HI Definition 2.1]. 



Definition 2.2 Let S safe/y (E0)-(E2). A weak ^-solution to fZ3] 
is a map 



it G C° (M + ;u + L 1 (M + ;R+ x IE 
«(t) G BV(R;i+xIR) for a.e. t £ 



(2.6) 



sitc/i i/iai 

(W) /or a// (/? G C*(IR + x M;M) whose support does not intersect x = 

(u c?t<^ + /(it) c^) dxdt = 0; 

for a.e. t G K + and u>ii/i \& as in ( fj?.3j) . £/ie junction condition is met: 
* (V,u(t,0-);a + ,u(i,0+)) =0. 

is a/so an entropy solution if 
(E) /or a// 9? G C£(IR + x R;IR + ) whose support does not intersect x = 

/ / (£(u) d t ip + F(u) d x tp) dxdt>0. 
ii+ii 

In the particular case of a Riemann Problem, i.e. of (jl.ip with initial datum 



u(0, x) 



n if x > 
u + if x < , 



Definition 12.21 reduces to [SJ Definition 2.1]. 

To state the uniqueness property in the theorems below, we need to 
introduce the following integral conditions, following [4, Theorem 9.2], see 
also [144 Theorem 8] and pQ. Given a function u = u(t,x) and a point (r, £), 
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we denote by Uf. T g, the solution of the homogeneous Riemann Problem 
consisting of (jl.2j> — fjl .3j) — A2.3H with initial datum at time r 

, . J lim x _^ u(r, x) if x < £ , . 

w(t,x) = < )(t a (2-7) 

lim. r __ > g + u(r, x) if x > 5 . 

and with £ satisfying (SO), (SI) and (S2). Moreover, define U) WT g\ as the 
solution of the linear hyperbolic Cauchy problem with constant coefficients 

I d t u + d x Au = t > t , 2 g s 

I w(t,x) = u(r,x) , 

with 2 = £>/ (n(r, £)) . 

The next theorem applies {$, Theorem 3.2] to ()1.2p with the choice (12. 3p 
to construct the semigroup generated by (|1.2|) - (|1.3|) - (|2.3p . The uniqueness 
part follows from [14} Theorem 2]. 

Theorem 2.3 Xei p satisfy (P) and S satisfy (^0)^(112). Choose any 
a > 0, u £ Aq. Then, there exist a positive A such that for all a~,a + with 
\a~ — a\ < A and \a + — a < A, i/iere exisi a map u as in \2. 5\) . positive 
5, L and a semigroup S: R + x T> —* 2? suc/i i/iai 

L P 5 {it £ u + L X (R; Aq): TV(u - u) < S}. 

2. For all u £ D, Squ = u and for all t,s > 0, StS s u = S s +t u - 

3. For all u,u' £ V and for all t, t' > 0, 

\\Stu — Sfu'\\ Ll < L ■ (\\u — u'\\ L i + \t — t'\ 

4- If u £ T> is piecewise constant, then for t small, Stu is the gluing of 
solutions to Riemann problems at the points of jump in u and at the 
junction at x = 0. 

5. For all u £ T>, the orbit t — > SfU is a weak ^-solution to |7j| ). 

6. Let A be an upper bound for the moduli of the characteristic speeds in 
B (u(R),<5). For all u £ V, the orbit u(t) = Stu satisfies the integral 
conditions 

(i) For allr > and £ £ R, 

I rt+hX 

lim-/ u(t + h, x) - Uf e Jr + h,x) dx = . (2.9) 
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(ii) There exists a C > such that for all r > 0, a,b G E and 
£ G ]a,6[, 



6- ft A 



a+/iA 



u(r + h, x) - U( U . T ^ (t + h, x) 



dx 

(2.10) 



< C 



TV{n(r);]a,6[} 



7. 7/ a Lipschitz map w: M — ► Z> satisfies $2.y\) - $2.1U\) . then it coincides 
with the semigroup orbit: w(t) = St (w(0)) . 

The proof is deferred to Paragraph 13. 11 Note that, similarly to what happens 
in the standard case of [U Theorem 9.2], condition (|2.10p is always satisfied 
at a junction. 

2.2 A Pipe with Piecewise Constant Section 

We consider now a tube with piecewise constant section 

71-1 

a = a X] — oo,xi\ + a nX[x n , 

3=1 

for a suitable n G N. The fluid in each pipe is modeled by (jl.2p . At each 
junction Xj, we require condition (II. 3D . namely 

for all j = 1, . . . , n, where 
y(a j - 1 ,uj;a j ,u+) = u ± = l im ^ Uj ( x ) , (2.11) 

We omit the formal definition of ^-solution to (jl.2h (jl.3)) in the present 
case, since it is an obvious iteration of Definition 12.21 



Theorem 2.4 Let p satisfy (P) and S satisfy (T I 0)-(T I 2). For any a > 

and any u G Aq there exist positive M,A,5,L,A4 such that for any profile 
satisfying 

(AO) a G PC (M; ]o - A, a + A[) with TV (a) < M, 
there exists a piecewise constant stationary solution 

71-1 

u = «0XJ-oo,«i[ + %X]^, %+1 [ + ti»X)*»,+oo[ 
i=i 

to lL2 ]) - \2.1i\) satisfying 

iij G TOf/i |uj — u| < <5 /or j = 0, . . . n 
^ (o 3 -_i ) u J -_i;Oj ) u J -) = for j = 1, . . . ,n 
TV(i)<MTV(fl) (2.12) 

and a semigroup S a :M + x 2? a — > 2? a smc/i i/iai 
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1. V a ^{u£u + L X (R; Aq): TV(u -u)<6}. 

2. S% is the identity and for all t,s>0, S?S% = S% +t . 

3. For all u, u' G V a and for all t, t' > 0, 



■i + \t 



4- If u G T> a is piecewise constant, then for t small, StU is the gluing of 
solutions to Riemann problems at the points of jump in u and at each 
junction Xj. 

5. For all u G V a , the orbit t — > Sfu is a weak ^/-solution to U.2\) - $2.1l\) . 

6. The semigroup satisfies the integral conditions 12.9\) - [2.10\) in 6. of 
Theorem [Q 



7. If a Lipschitz map w: R — > T> satisfies \2.9\) - [2.10\) . then it coincides 
with the semigroup orbit: w(t) = St (w(0)) . 

Remark that 5 and L depend on a only through a and TV(a). In particular, 
all the construction above is independent from the number of points of jump 
in a. For every u, we provide below an estimate of M at the leading order 
in 5 and A, see p. lip and (|3.8[) . In the case of E as in (|2.19p and with the 
isothermal pressure law, which obviously satisfies (P), 



P(p) 



c 2 p. 



(2.13) 



the bounds (13. lip and (|3.8p reduce to the simpler estimate 



M 



a 

Ae~ 

a I - (v/cf 
4e (v/c) 2 



ifv/c G 0,l/\/2 
i£v/c G ll/V2,l 



(2.14) 



where v = q/p. Note that, as it is physically reasonable, M is a weakly 
decreasing function of v, so that at lower fluid speeds, higher values for the 
total variation of the pipe's section can be accepted. 

Furthermore, the estimates proved in Section 13.21 show that the total 
variation of the solution to (|1.2|) - (|2.11|) may grow unboundedly if TV (a) is 
large. Consider the case in Figure HJ A wave cr^~ hits a junction where the 
pipe's section increases by Ao > 0. From this interaction, the wave crj of 
the second family arises, which hits the second junction where the section 
diminishes by Aa. At the leading term in Aa, we have the estimate 

< I l + K(v/c) ' 



(To 




f7n 



where 



2(l-£) 3 (l + £) s 



(2.15) 
(2.16) 
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a 



Aa 




Figure 1: A wave a 2 hits a junction, giving rise to a 2 which hits a second 
junction. 



see Section [3721 for the proof. Note that /C(0) = —1 whereas lim^_>i_ /C(£) = 
+oo. Therefore, for any fixed Aa, if v is sufficiently near to c, repeating 
the interactions in Figure [TJ a sufficient number of times makes the 2 shock 
waves arbitrarily large. 

2.3 A Pipe with a W 1 1 Section 

In this paragraph, the pipe's section a is assumed to satisfy 

{a e W 1 ' 1 (R; ]a — A,a + A[) for suitable A > 0, a > A 
TV(o) < M for a suitable M > 
a'(x) = for a.e. i£l\[-I,I] for a suitable X > . 

For smooth solutions, the equivalence of (jl.lj) and (jl.4[) is immediate. Note 
that the latter is in the standard form of a ID conservation law and the 
usual definition of weak entropy solution applies, see for instance |19[ Defi- 
nition 3.5.1] or [9] Section 6]. The definition below of weak entropy solution 
to (|l.ip makes the two systems fully equivalent also for non smooth solutions. 

Definition 2.5 A weak solution to (jl.ip is a map 
u £ C° (r + ; u + L^R; R + xI 
such that for all v? G Cj(M + x R;R) 



ap 
aq 



dtp + 



aq 
aP(u) 



d x sp + 





p(p)d x a 



(p\dxdt = 0. (2.17) 
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u is an entropy weak solution if, for any tp G C£(IR + xR;l), <p > 



JM.+ JR 




(2.18) 



Lemma 2.6 Let a satisfy (Al). Then, u is a weak entropy solution to U.l\) 
in the sense of Definition \2. 51 if and only if it is a weak entropy solution 



The proof is deferred to Section 13.31 

Now, the section a of the pipe is sufficiently regular to select stationary 
solutions as solutions to either of the systems fjl .5j) , which are equivalent by 
Lemma 12.61 Hence, the smoothness of a also singles out a specific choice of 
S, see [3 formula (14)]. 

Proposition 2.7 Fix a~,a + £ ]o — A,o + A[ and u~ £ Aq. Choose a 
function a strictly monotone, in C 1 , that satisfies (Al) with a(—X—) = a~ 
and a{X+) = a + . Call p = R a {x;u~) the p-component of the correspond- 
ing solution to either of the Cauchy problems lll.5\) with initial condition 
u{—X) = u~ . Then, 

1. the function 



satisfies (E0)-(E2); 

2. if a is a strictly monotone function satisfying the same requirements 
above for a, the corresponding map E coincides with E. 

The basic well posedness theorem in the present W ' case is stated similarly 
to Theorem 12.41 

Theorem 2.8 Let p satisfy (P). For any a > and any u G Aq there exist 
positive M,A,d,L such that for any profile a satisfying (Al) there exists a 
stationary solution u to satisfying 




(2.19) 



u £ Aq with \\u(x) — u\\ < 5 for all x Gl 



and a semigroup 5 a :M + x T>' 



<a 



V a such that 



1. V a ^ {ueu + L X (R; Aq): TV(u — u) <6}. 




is the identity and for all t, s 



> 0, s?s* = si 
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3. for all u, v! G V a and for all t, t' > 0, 

\\S?u - S?,u'\\ L1 < L • (\\u - u'\\ L1 + \t- t'\\ . 

4- for all u E V a , the orbit t — > Sfu is solution to in the sense of 

Definition \2.5[ 



5. Let A be an upper bound for the moduli of the characteristic speeds in 
B (-u(R),(5). For all u G T>, the orbit u(t) = StU satisfies the integral 
conditions 

(i) For allr < and £ G R, 



lim — 

h->o h 



£-hX 



u(t + h,x) -UL t ^Jt + h,x) dx = 0. (2.20) 



(ii) There exists a C > such that, for all r > 0, a,b G R and 

Z e KM, 



b-hX 



a+hX 



u(t + h,x)- U (u . T ^ (t + h, x) 



dx 



(2.21) 



< C 



TV{«(r);]a,6[}+TV{a;]a,6[} 



6. If a Lipschitz map w: R — > P solves lll.l]) , then it coincides with the 
semigroup orbit: w(t) = St (w(0)) . 

Thanks to Theorem E3J the proof is obtained approximating a with a piece- 
wise constant function a n . The corresponding problems (|1.2p - (|2.1ip gener- 
ate semigroups defined on domains characterized by uniform bounds on the 
total variation and with a uniformly bounded Lipschitz constants for their 
time dependence. Then, we pass to the limit (see Section 4 for the proof) 
and we follow the same procedure as in [3J Theorem 9.2] and theorems 2 
and 8] to characterize the solution. 

As a byproduct of the proof of Theorem 12.81 we also obtain the follow- 
ing convergence result, relating the construction in Theorem 12.41 to that of 
Theorem [~~ 



Proposition 2.9 Under the same assumptions of Theorem \2.8l for every 
n E N, choose a function j3 n such that: 

(i) P n is piecewise constant with points of jump y^, ...,y™ n , with y\ = 
-X, y™« = X, and max^j/^ 1 - y J n ) < l/n. 

(ii) p n (x) = for all x G R\ [-X,X]. 

(Hi) P n — > a' in L 1 (R;R) with H^Hl 1 — M, with M as in Theorem [2~. 
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Define a n {x) = a(— X— ) + JZ x f3 n (£)dt; and points x 3 n G 
j = 1, . . . , m n — 1 and let 

rn n - 1 

a n = a(-X-) X] + ^ an (^+i) x ^^ +1[ + a (X+) X[ ^ i+oo[ 

i=l 

(see Figure^) . Then, a n satisfies (AO) and the corresponding semigroup 
S n constructed in Theorem \2.J\ converges pointwise to the semigroup S con- 
structed in Theorem \2.8l 



j i+1 



for 



3 Technical Proofs 

3.1 Proofs Related to Section D 

The following equalities will be of use below: 

-Ai A2 and 



d p P 



d q P = Ai + A 2 . 



(3.1) 



Proof of Lemma \2.1i Apply the Implicit Function Theorem to the equal- 
ity f = in a neighborhood of (a,u,a,u), which satisfies ^ = by (SI). 
Observe that d u T,(a, a;u~) = by (SI). Using (|3,1|) . compute 



det d u +^(a, u, a, u) 



det 



det 



— dp+T.1 a 

a d p + P a d q + P 

a 

a d p + P a d q + P 



= a 2 Xi(u) \ 2 (u) 
+ 0, 



completing the proof. 



□ 



Proof of Theorem \2.3l Let A be defined as in Lemma 12.11 Assump- 
tion (F) in Theorem 3.2] follows from (P), thanks to (l2~Tj) and to the 
choices f|2 . 2 [) — f|2 . 5 [) . We now verify condition [8j formula (2.2)]. Recall that 
-D u -£(a, a; u) = by (SI). Hence, using (|3.ip . 



det [D u -^(a, u; a, u) ■ r\(u) D u +ty(a, u; a, u) • r 2 (u)~\ 

aAi(u) + d p -T,i(a, a; u) + X\d q - Si (a, a; u) aX 2 (u) 

a (Ai(u)) 2 + <9 p -£ 2 (a, a;u) + X^d q -T, 2 (a, a; u) a (X 2 (u)) 2 

a\i(u) aX 2 {u) 



det 



a (Xi(u)) z a (X 2 (u)Y 
= a 2 Xi(u) X 2 (u) (X 2 (u) - X\(u)) 
+ 0. 
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The proof of 1.-5. is completed applying [8, Theorem 3.2]. The obtained 
semigroup coincides with that constructed in Theorem 2], where the 
uniqueness conditions 6. and 7. are proved. □ 



3.2 Proofs Related to Section D 

We now work towards the proof of Theorem 12.41 We first use the wave 
front tracking technique to construct approximate solutions to the Cauchy 
problem (|l,2p - (|2.11|) adapting the wave front tracking technique introduced 
in [U Chapter 7]. 

Fix an initial datum u Q E u + L^IRj^o) and an e > 0. Approximate 
u with a piecewise constant initial datum u £ Q having a finite number of 
discontinuities and so that lim e ^o \ \ u o ~ n o|| L i =0. Then, at each junction 
and at each point of jump in u E D along the pipe, we solve the corresponding 
Riemann Problem according to Definition 12.21 If the total variation of the 
initial datum is sufficiently small, then Theorem 12.31 ensures the existence 
and uniqueness of solutions to each Riemann Problem. We approximate 
each rarefaction wave with a rarefaction fan, i.e. by means of (non entropic) 
shock waves traveling at the characteristic speed of the state to the right of 
the shock and with size at most e. 

This construction can be extended up to the first time t\ at which two 
waves interact in a pipe or a wave hits the junction. At time i\ the functions 
so constructed are piecewise constant with a finite number of discontinuities. 
At any subsequent interaction or collision with the junction, we repeat the 
previous construction with the following provisions: 

1. no more than 2 waves interact at the same point or at the junction; 

2. a rarefaction fan of the i-th family produced by the interaction between 
an i-th rarefaction and any other wave, is not split any further; 

3. when the product of the strengths of two interacting waves falls below 
a threshold e, then we let the waves cross each other, their size being 
unaltered, and introduce a non physical wave with speed A, with A > 
sup( M ) A2 (it); see |U Chapter 7] and the refinement [2]. 

We complete the above algorithm stating how Riemann Problems at the 
junctions are solved. We use the same rules as in § 4.2] and [HJ § 5]. In 
particular, at time t = and whenever a physical wave with size greater 
than e hits the junction, the accurate solver is used, i.e. the exact solution 
is approximated replacing rarefaction waves with rarefaction fans. When a 
non physical wave hits the junction, then we let it be refracted into a non 
physical wave with the same speed A and no other wave is produced. 

Repeating recursively this procedure, we construct a wave front tracking 
sequence of approximate solutions u £ in the sense of [H Definition 7.1]. 
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At interactions of waves in a pipe, we have the following classical result. 



Lemma 3.1 Consider interactions in a pipe. Then, there exists a positive 
K with the properties: 

1. An interaction between the wave of the first family and a 2 of the 
second family produces the waves at and a 2 with 



+ 



ff, — (To 



< K ■ 



°1 °2 



(3.2) 



2. An interaction between a[ and a'( both of the same i-th family produces 
waves of total size af and a 2 with 



+ 



o- 



+ 



°2 ~ (°2 + 02) 



< K ■ lo-X'l 



<K-\a' 2 a'i\ 



ifi = l, 
ifi = 2. 



3. An interaction between the physical waves o~ 1 and a 2 produces a non 
physical wave , then 



(To 



< K ■ 



4- An interaction between a physical wave a and a non physical wave o" 3 
produces a physical wave a and a non physical wave a^, then 



(To 



fTo 



< K ■ 



For a proof of this result see [H Chapter 7] . Differently from the construc- 
tions in [3 [8] , we now can not avoid the interaction of non physical waves 
with junctions. Moreover, the estimates found therein do not allow to pass 
to the limit n — > +oo, n being the number of junctions. 

Lemma 3.2 Consider interactions at the junction sited at Xj. There exist 
positive K\, K 2 , K% with the following properties. 



1. The wave <r 2 hits the junction. The resulting waves o~^,o~ 2 satisfy 







\ a j~ a J-l\ 


a 2 



0~n 



< 1 1 + K2 \aj — cij-i 



(To 



< e 



K2 \ CLj —0, 



3 a 3 — 1 



(To 
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2. The non-physical wave a hits the junction. The resulting wave o" H 
satisfies 



a 



< (1 + K 3 \a,j - dj-i 

< e K A a i~ a i-A 



a 



u 



a 



(7 



a 



Proof. Use the notation in the figure above. Recall that af and a 2 
are computed through the Implicit Function Theorem applied to a suit- 
able combination of the Lax curves of (jl.2p . see Proposition 2.4] and [5J 
Proposition 2.2]. Repeating the proof of Theorem 12.31 one shows that the 
Implicit Function Theorem can be applied. Therefore, the regularity of the 
Lax curves and (P) ensure that crj*" = \ (J 2^ a i ~ a j-i'^) an d a t = 



0~ o Cn , Cti 



). An application of [H Lemma 2.5], yields 



af (0, aj — dj-i) uj 



< Ki \a 



i\aj — aj-i\ 



a,, 



at (0, a,- — a-i 



'2 V u ' "J ~~ u 'j-l) u ) 



Or, - (Jo 



(T 



< 



< Qj — a i-i| 



1 + K 2 \aj — a,_i 



(To 



completing the proof of 1. The estimate at 2. is proved similarly. 



(3.3) 



□ 



We now aim at an improvement of (j3.3f) . Solving the Riemann problem 
at the interaction in case 1. amounts to solve the system 



(3.4) 



£ 2 T [Ci{u;a{) \ ;a : 



T [C 2 (u;a 2 



By (|2.ip . the first order expansions in the wave's sizes of the Lax curves 
exiting u are 



C\{u;a) 



p - a + o(a) 
q - Xi(u)a + o(a) 



and jC,2(u;ct) 



p + a + o(a) 
q + A 2 (n)cr + o(a) 



while the first order expansion in the size's difference Ao = a + — a of the 
map T defined at (|2.4f) . with v = q/p, is 



T(a,a + Aa;w) 



if 



l + #^jp + (Aa) 
1 - 4^ <7 + o(Ao) 



where (3.5) 



2 , d a+ E-p(p) 
C 2 — V 2 
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Inserting these expansions in (13, 4D . we get the following linear system for 



i 



Aa 



(7 



1 + i?4r ) 



: C 



where 



if 



v 2 + (9 a +S(a,a,u) -p(p)) /p 



c 2 — u 2 



and G 



(c'(p)p — v)H — v 
V + c 



and all functions are computed in tt. The solution is 



0V 



(To 



.* ( i +fl+ *)£,; 

' x _ Ai^ + A 2 (l + G) Aa 
2c a 



(3.6) 
(3.7) 



which implies the following first order estimate for the coefficients in the 
interaction estimates of Lemma 13.21 



Ki = — 
2a 



K, 



1 

2a 



<9 Q+ E-p( P ) 

p 



l-2(f) 2 + ^(^) 2 + ^(^-l 



9 a+ E-p(p) 



(3.8) 



The estimate (|3,7p directly implies the following corollary. 

Corollary 3.3 If aj — 2s sufficiently small, then and are ei- 

ther both rarefactions or both shocks. 

Denote by a\ the wave belonging to the i-th family and sited at the 
point of jump x a , with x a in the j-th pipe Ij, where we set Iq = ] — oo,xi[, 
Ij = ]xj, Xj + \ [ for j = 1, . . . , n — 1 and I n = ]x n , +oo[. Aiming at a bound 
on the Total Variation of the approximate solution, we define the Glimm-like 
functional, see [S formulae (7.53) and (7.54)] or also [El 13 QU |20] , 



V = E E (\cri a e c ^\ a »- a »-i\ + 

j=0 x a £lj 



'2.a 



CY7h=\ | a h+l — a h\ 



j=0 



cZZZ-\*h + i-a h \ \a\, 
a non physical in Ij 
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Q 



E 

V + Q: 



3 3 
a i,a a i>,a' 



(3.9) 



where C is a positive constant to be specified below. A is the set of pairs 
(°i a ' °f' «' ) °f approaching waves, see [3, Paragraph 3, Section 7.3]. The 
i-wave a\ a sited at x a and the i'-wave o\, , sited at x Q ( are approaching 

if either i < i' and x a > x a >, or if i = i' < 3 and min{<7£ a > cr|, ,} < 0, 
independently from j and j'. As usual, non physical waves are considered 
as belonging to a fictitious linearly degenerate 3rd family, hence they are 
approaching to all physical waves to their right. 

It is immediate to note that the weights exp [cY^h=i \ a h ~ a h-i\\ an d 
,; j ~ a h\) in the definition of V are uniformly bounded: 



exp (CY2,=) K+i 



1 < exp [CJ2i=i \ a h ~ o/i-il 



1 < exp(CE n_1 



< exp(CTV(a)) 

< exp(CTV(a)) 



(3-10) 

Below, the following elementary inequality is of use: if a < b, then 
e a - e b < -(b - a)e a . 



Lemma 3.4 There exists a positive 5 such that if an e- approximate wave 
front tracking solution u = u(t, x) has been defined up to time i, T (u(i— )) < 
5 and an interaction takes place at time i, then the e-solution can be extended 
beyond time i and T (Wi+)) < T (u(i—)). 



Proof. Thanks to (I3,10p and Lemma 13.11 the standard interaction esti- 
mates, see [U Lemma 7.2], ensure that T decreases at any interaction taking 
place in the interior of Ij, for any j = 0, . . . , n. 

Consider now an interaction at Xj. In the case of[T]in Lemma l3.2( 



< 



AQ 



E 



0~\ Oi, a 



+ 



E 

(crJ,tTi, a )€.4 



< Ki \<lj - Oj_i| ^ \o-j ta \ + 



< (Ki + K 2 ) T(t-) \ aj - Oj-i | 



< (K\ + K 2 ) 5 \dj — dj-i | 

Ay 



(T 



3 Jf 2 Oj— Oj_l I 



t.Q 
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< e C T, J hJi\ a h-a h -i\ 



_|_ gCELj | a h+i- a /i| 



CTn 



CYX=}-1 \ a h+l-1h\ 



(J 



— O7 — 1 



CELi | a h+l~ a h| g-^2|oj— _ e CJ2h=j-l \ a h+l— 0-h 



a,, 



CTn 



_|_ gCEftlj |lh+l-1h| ^glf2[Oj--Oj_i| _ e C|ttj— Oj_i 



-(C - Er 2 )| aj - 0i _i| e^l^'-^-il e CE ^ l a *+i-«fc| 



CTn 



< 



< 







\ a i ~ a j-i\ 


a 2 



AT 



(iTi + if 2 ) (l + e ^l° + -° I) e CTV M + (iTi + K 2 )5 - C^j 



Choosing now, for instance, 
1 



S< 1, C 



TV(a) 



a — a 



< |? and TV(.) < I (3.11) 



the monotonicity of T in this first case is proved. 

Consider an interaction as in[2j of Lemma 13.21 Then, similarly, 



AQ < 

0+ ,<Ti, a )eA 







I "*i,Q! | 


(T 



< T(t— ) |aj — a,j-\\ |cr 



(7 



(7 



< ^£3 8 \aj — dj-x\ 
AV < e c ^h^\ ah + 1 ~ ah 

< ^ C ^h=j\ a h+l-<*h\ e K3\a.j-aj-i\ _ & C Eh=]_i \ a h+i -a h 

< e c TZZ)\ a h+i-a h \ f e K 3 \a j -aj-i\_ e C\aj-a j -i 

< (K 3 - C)\ aj - Oj_i| e^h^-il e CE ^ 
AT < (^ 3 e X3 l a+ - a "le CTV W+^-C) | % - - a,_ a 



a 

-an 



a 



a 
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and the choice 5 < 1 and C > ensures that AT < 0. 



□ 



Proof of Theorem \2.4[ First, observe that the construction of the station- 
ary solution u directly follows from an iterated application of Lemma 12.11 
The bound (|2.12p follows from the Lipschitz continuity of the map T defined 
in Lemma 12.11 Define 

T> = |u G u + 1 L 1 (M.;A ):u G PC and T(u) < , 

where PC denotes the set of piecewise constant functions with finitely many 
jumps. It is immediate to prove that there exists a suitable C\ > such that 
^TV(u){t) < V(t) < CiTV(u)(i,-) for all (it) G V. Any initial data in V 
yields an approximate solution to (ll.2p attaining values in T> by Lemma [3. 41 
We pass now to the L 1 -Lipschitz continuous dependence of the approx- 
imate solutions from the initial datum. Consider two wave front tracking 
approximate solutions u\ and u 2 and define the functional 



EE 

3=1 i=l ' 



r+oo 






s{{x) 


lo 





(3.12) 



where s](x) measures the strengths of the i-th shock wave in the j-th pipe 
at point x (see [U Chapter 8]) and the weights W- are defined by 



wiix) 



1 + ki A{(x) + m k 2 (T(ui) + T(u 2 )) 



for suitable positive constants k±, k 2 chosen as in [;4, formula (8.7)]. Here T 
is the functional defined in (|3.9p . while the A\ are defined by 



A j A 




x a < x, a G Jj(u\) 
x a > x, a G Jj{u 2 ) 

x a < x, a G Jj{u 2 ) 
x a > x, a G Jj{u\) 



if (x) < 0, 



if sf(x) > 0; 



see [U Chapter 8]. Here, as above, a\ a is the wave belonging to the i-th 

family, sited at x a , with x a G Ij. For fixed ki, k 2 the weights VF/(x) are 
uniformly bounded. Hence the functional $ is equivalent to L 1 distance: 



1 

C~2 



\ui - n 2 || L i < $ {ui,u 2 ) < C 2 • - n 2 || L i 
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for a positive constant C2. The same calculations as in [U Chapter 8] show 
that, at any time t > when an interaction happens neither in u\ or in U2, 



dt 



*(t*i(t),U2(t)) < C 3 e 



where C3 is a suitable positive constant depending only on a bound on the 
total variation of the initial data. 

If t > is an interaction time for u\ or U2, then, by Lemma 13.41 
A T + T (u2(t)) < and, choosing K2 large enough, we obtain 

A$ («i(i),«2(t)) <0. 

Thus, $ (ui(i),n 2 (t)) - $ («i(s), u 2 («)) < C 2 e (t - s) for every < s < t. 
The proof is now completed using the standard arguments in [H Chapter 8] . 

The proof that in the limit e — > the semigroup trajectory does indeed 
yield a ^-solution to (|1.2p and, in particular, that (|2.11|) is satisfied on the 
traces, is exactly as that of [U Proposition 5.3], completing the proof of 1.-5. 

Due to the local nature of the conditions (|2.9p - (|2.10p and to the finite 
speed of propagation of (|1.2p . the uniqueness conditions 6. and 7. are proved 
exactly as in Theorem 12.31 □ 



Proof of estimate ^2.1J$ . We first compute d a +H, with S defined 
in (|2,19p . To this aim, by 2. in Proposition 12.71 (in Paragraph 2.3), we 
may choose 



a(x) 



if x € ]— 00, —X[ , 
-(x + X) + cT if x 6 [~X,X], 



a + if x S ]X, +00 [, 

so that we may change variable in the integral in (|2.19|) to obtain 

d a+ ^ = d a+ ^ p{R a (a,u)) daj = P (p) + 0(Aa). 



(3.13) 



Now, estimate (12T41 directly follows inserting (12TT3I) and (13TT31 in (|3TLj) 
and ([SID. □ 

Proof of estimates ^2.15\) - $2.16\) . Refer to the notation in Figure [U 
where the pipe's section is given by 



a\x) 



a if x G ]— 00, Z[, 

a + Aa if x G [Z, 2/] , 
a if x G ]2/,+oo[, 
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where Aa > 0. The wave a 2 arises from the interaction with the first 
junction and hence satisfies (|3.7|) . Using the pressure law (|2.13p and (|3,13p . 
we obtain 



a 2 = (l + ip (u, a) Aa) a 2 , where 



1 



1/2 



1 - (v/c) 2 



Now we iterate the previous bound to estimate the wave a 2 + which arises 
from the interaction with the second junction, i.e. 



ovJ~ + = (l - ip(a + Aa, u + ) Aaj a£, 



where, by (|3.5p . 

ifi(a + Aa, u + ) = ifi 



^ . ( 1 Aa\ / Aa\ ^ 
a + Aa, I 1 + — — | p, ( 1 - — ) q 



V 



/ 



Introduce rj = 1/ (l — (v/c) 2 ) and 7? = Aa/a to get the estimate 



a. 



++ 



1 + (if}(a, u, ) — ip(a + Aa, u + )^ Aa^j a 2 



Aa ( r]\ Aa 

1 + ( — 1 + 77 ] + t- 

a V 2 J a + Aa 



1/2 



V 

/ / 

Aa 



1-1? V 



<j 



n 



i 



V 



77 1 
+ X + 



2 1 + 7? 



V 



1/2 



1-1? u 



cr 



and a further expansion to the leading term in Aa gives (|2.15p -( T2.16p . □ 



3.3 Proofs Related to Section 

Proof of Lemma XZM. If a G C 1 (R; [a~,a+]) and u is a weak entropy 
solution of (11.4|) . Then, 



9 



ap x 



(^1 (it 
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/ / 

/R+ JR \ 


ap 
aq 


a 


aq 
aP(u) 


/ /f 

m+ Jr \ 


ap 
aq 


d t * + 
a 


aq 
aP(u) 



-d x (f 



d x * + 
a 



aq 






■) 


p(p)d x a 





d T a dx dt 



dx dt 



showing that (|2.17p holds. Concerning the entropy inequality, compute pre- 
liminarily 



V(aE{u)) 



a x 

^d x a 

ap x 



2p 2 
2p< 



+ 



+ 



ld r a 

y 

ap '" 



^-d T a 



p * P( r ) , . Q i s , T \ o 
—^dr + -p(p) + -^)d x a 



= ±(E(u)+p(p)) 8 x a 
p 

= F(u) d x a . 

Consider now the entropy condition for (jl.4p and, by the above equality, 



< / / ( E{u) d t <p + F(u) d x tp - VE(u 
Jr+Jr \ 

(aE(u) d t - + aF(u) d x - 
. \ a a 

+ F(u)d x a - V(aE(u)) 



\d x a 
^d x a 

ap x 



if I dxdt 



7;d x a 
ir,d x a 

ap x 



dx dt 



[ [ (aE{u)d t ^ + aF{u)d x ^ 

JR+ JR K a a 



+ (F(u)d x a - F(u)d x a) -)dxdt 

' a) 



= [ [ (aE(u)d t - + aF(u)d x -) dxdt, 
Jr+ Jr V a a / 

showing that (I2,18p holds. The extension to a 6 W 1 ' 1 is immediate. 



□ 



Proof of Proposition \2.1\ The regularity condition (£0) follows from 
the theory of ordinary differential equations. Condition (SI) is immediate. 

Consider now the item 2. If a\ and 02 both satisfy (Al), are strictly 
monotone, smooth and have the same range, then a\ = a2°(f for a suitable 
strictly monotone (p with, say (p' > 0, the case (p' < is entirely similar. 
Note that if u = (Ri(x;u~),Qi(x;u~)) solves (II. 5p with a = at, then di- 
rect computations show that Ri(x,u~) = R2 (tp(x),u~^ and Qi(x,u~) = 
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Q2(ip(x),u ). Hence 



x 



Si (a , a + ;u ) = J p^R\(x;u )j a'i(x)dx 

= J V \R2 (p(x);u~yj a' 2 (<p(x)) (p'(x)dx 

= S 2 (a~,a + ;ti _ ) . 

Having proved (SO) and (SI), we use the map T defined in Lemma 12.11 
We first prove that S satisfies £(a~, a + ; u~) + £ (a + , a - ; T(a + , a~; u~)) = 0, 
given a satisfying (Al), strictly monotone and with a(—X) = a~ , a(X) = 
a + , let a[x) = a~ + a + — a(x). Then, using 2. proved above, and integrat- 
ing (II. 5p backwards, we have 

^ (^a + ,a~ ]T{a~ ,a + ■ 1 u~)^ = J p (^R{x]T{a~ ,a + ]u~)^ a(x)dx 

p (^R(x; a~ , a + ; a'{x)dx 



X 
X 



= — S(a , a ; u ) . 

Finally, condition (S2) follows from the the flow property of R and the 
additivity of the integral. Indeed, by 2. and 3. we may assume without loss 
of generality that a~ < a < a + . Then, let q = Q(x; u~) be the q component 
in the solution to (jl.5p with initial condition u(0) = u~. Then, if T is the 
map defined in Lemma 12. 11 we have 

T(a~,a + ;u~) = (^i?(a" 1 (a + ); if"), Q(a _1 (a + ); u~)) 

so that 

S(a _ , a + ; u~) 

p^i?(x,u~)^ a'(x)dx 

r~ 1(a0) / -\ / ■ _ \ 

/ p I i?(x, « ) a (x) dx + / p I i?(x, « ) a (x) 



A' 



a" 1 (a ) 



p(^R(x,u )J a!{x)dx 
f p(R(x,R(a- 1 (a ),u-),Q(a- 1 (a°),u^))^) a'{x)dx 



l(a°) 

S(a _ , a ; u~) + S ( a , a + ; T(a~, a + ; u~) 
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proving 1. 



□ 



Proof of Theorem [Ql Fix a > 0, and u £ A . Choose M,A,L,S as 
in Theorem 12.41 With reference to these quantities, let a satisfy (Al). For 
n G N, let a n ,a n ,[3 n be as in Proposition 12.91 Note that a n is piecewise 




Figure 2: The thick line is the graph of a = a(x), the dotted line represents 
a n while the polygonal line is a n 

linear and continuous. By (iii), we have that a n — > a and a n — > a in L 1 . 
Moreover, TV(a n ) < M and TV(a n ) < M and, for n sufficiently large, 
flnO^) ^ ]a — A,a + A[. Hence, for n large, a n satisfies (AO). Call S n the 
semigroup constructed in Theorem 12.41 and denote by T> n its domain. 

Let be a sequence of initial data in T> n . The S n are uniformly Lipschitz 
in time and S^u^ have total variation in x uniformly bounded in t. Hence, 
by [H Theorem 2.4], a subsequence of u n {t) = Sfu^ converges pointwise 
a.e. to a limit, say, u. For any ip G C*(R + x and for any fixed n, let 

e > be sufficiently small and introduce a C£°(R;]R) function rj e such that 



T] £ {X) 
T) £ (x) 

Thus, we have 



for all x G {jj={ l [x{ - e,x{ + e] , 

1 for all x G \J^~ 2 [xi + 2e,x j n +1 - 2e\ 








( 


O'nPn 




Q"nQn 


( 


O'nPn 




QnQn 



d x (p dx dt 



Ve d t (p + 



dnQn 

a n P{u n ) 

QnQn 

a n P(u n ) 



r] £ d x ip dx dt 



d x (Ve <p) dxdt 
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lim 





O'nQn 


/ / 

m+Jm 


a n P(u n ) 



(fd x rj £ dx dt . 



The first summand in the latter term above vanishes by Definition 12.21 ap- 
plied in a neighborhood of each x- n . The second summand, by the BV 
regularity of u n , converges as follows: 



lim 

m n — 1 

E 

3=1 
m n — l 

E 

i=i 




O'nPn 

a n P(u n 



O'nQ.n 

a n P(u r , 



d t <f + 

ip d x rj £ dx dt 



d x (p dx dt 



On(^'n~t _ )<J l n(2'n~l _ ) S n (li ^)Qn{%n ) 
Q'n(%n~\~)Pn(.%n~\~) Cbn^n )Pn(%n ) 



E I d n {Xn~ ), CL n {Xn~\~) , u(t, Xn~) 



ip(t,x J n )dt 
ip{t,x J n )dt. 



We proceed now considering only the second component. Using the map 

mn — l 

<p n (t,x) = <p(t,x)x ] _ 0Otyk[ (x)+ J2 V&^X^yi+lfe) 



i=i 



+<p(t,x) X]y ^ )+oo[ (x), 



we obtain 



m n — 1 

E 

i=i 

m n — 1 

E 

3=1 
rrin—1 

E 

3=1 
mn—l 



E 



S \ a n (x J n -),a n (x J n +),u(t, x 3 n -)j <p(t, x J n ) dt 
s [a n (yi),a n (yi +1 ),u(t,x 3 n -)j (p(t,x J n )dt 

s (anM),an(l/i + VMn-)) vfax-^dt 

p( R an (^x; u n (t, x^—)^ j a' n {x) dx ip(t, x 3 n ) dt 

1 / \ 
p[ R an (x; u n (t, x J n —) ) I a' n (x) dx (p(t, x J n ) dt 



v j+1 



3=1 Vr 
m.n — 1 

p( Ra " [x;u n (t,xi-) 

[ 3=1 
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xa' n {x)tp{t,x 3 n )x [y i it yj+i [ 



(x) dx dt 




p(p{xj) d x a{x) (p{t,x 3 n ) dx dt 



as n 



+00 



where we used (i) in the choice of the approximation a n . 

We thus constructed a solution to (jl.ip . for any initial datum in T>. Note 
that this solution satisfies (I2,20p - (|2.2ip . as can be proved using exactly the 
techniques in [14, Theorem 8]. Therefore, the whole sequence u n converges 
to a unique limit u, which is Lipschitz with respect to time. This uniqueness 
implies the semigroup property 2. in Theorem 12.81 The Lipschitz continuity 
with respect to the initial datum follows from the uniform Lipschitz regu- 
larity of the approximate solutions u n , completing the proof of 3. Finally, 
6. is proved exactly as in [14; Theorem 8]. □ 
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